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Quark counting rules derived from recombination models agree well with data on hadron produc¬ 
tion at intermediate transverse momenta in relativistic heavy-ion collisions. They convey a simple 
picture of hadrons consisting only of valence quarks. We discuss the inclusion of higher Fock states 
that add sea quarks and gluons to the hadron structure. We show that, when recombination oc¬ 
curs from a thermal medium, hadron spectra remain unaffected by the inclusion of higher Fock 
states. However, the quark number scaling for elliptic flow is somewhat affected. We discuss the 
implications for our understanding of data from the Relativistic Heavy Ion Collider. 


The production of hadrons with transverse momenta 
of a few GeV/c in Au-I-Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) has been found to exhibit 
certain unusual features. The emission of baryons, in 
comparison with mesons, is significantly enhanced. For 
example, the ratio of protons to pions is about three times 
larger than in p -|- p collisions [ij. Furthermore, the nu¬ 
clear suppression factor Raa below 4 GeV/c is close to 
unity for protons and A-hyperons, while it is about 0.3 
for pions The azimuthal anisotropy (called “elliptic 
flow” V 2 ) of baryons lags behind that of mesons at low 
transverse momenta but exceeds the anisotropy of me¬ 
son emission above 2 GeV/c 0, 0. The recombination 
of quasi-thermal, deconfined quarks has been proposed as 
an explanation for these peculiarities 

These recombination models are based on the concept 
of constituent quark recombination, which assumes that 
the probability for the emission of a hadron from a de¬ 
confined medium is proportional to the probability for 
finding the valence quarks of the hadron in the density 
matrix describing the source. The baryon enhancement, 
as well as the different momentum dependence of meson 
and baryon anisotropies, rely essentially on the different 
number of valence quarks in mesons (two) and baryons 
(three). The simplicity of this concept has been criti¬ 
cized, because it does not do justice to the complexity 


of the internal structure of hadrons in quantum chromo¬ 
dynamics (QCD). It is the purpose of this article to ad¬ 
dress some of these objections by showing that important 
features of the emission of relativistic hadrons by sud¬ 
den recombination from a thermal deconfined medium 
do not depend on specific assumptions about their inter¬ 
nal structure and survive in a more complete description 
of the emitted hadrons. 

The quark-gluon content of a hadron in QGD, as that 
of a bound state in any strongly interacting quantum 
field theory, is a function of the resolution scale . The 
larger the scale is, the more quantum fluctuations can 
be seen, and the number of gluons and sea quarks in a 
hadron increases. Moreover, the structure of a hadron 
depends on the frame of reference. This is so because the 
boost operator K and the Hamiltonian H generally do 
not commute: [K, H] ^ 0. Furthermore the quark-gluon 
structure depends on the measurement process, i. e. it is 
not universal. Different processes probe the Fock states 
in a hadron with different weights. 

In the following we will use a light-cone frame, where 
formally the hadron momentum P —> 00 and the mo¬ 
mentum fractions of the partons are the only dynamic 
degrees of freedom. A meson M with valence quarks qa 
and qfs can then be written as an expansion in terms of 
increasingly complex Fock states: 


\M) = / dXadXbS{Xa+Xb-l)ci{Xa,Xb)\qaiXa)qi3iXb)) 


dXadXbdXcS{Xa + Xb+Xc- l)c2{Xa, Xb, Xc) \qa{Xa)qi3{xb)giXc)) 


( 1 ) 


dxadxbdxcdxdS{xa + Xb + Xc + Xd - l)c3(a;a, Xb, Xc, Xd) \qa{xa)q0{xb)q-y{xc)q-y{xd)) + 


Hence XiP is the momentum of parton i. We have sup- pressed the momentum P in the notation for simplicity. 
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For convenience, we adopt a probabilistic normalization 
for partons: {qa{xa)\qi3{.Xb)) = 5ap5{xa - Xb). This im- 
plies that the completeness relation for the Fock space 
expansion takes the simple form: 


1 — Cl + C 2 + C 3 


= / dXadXb5{Xa + Xb - 1) \Ci{Xa,Xb)\" 

Jo 


( 2 ) 


+ 

+ 



\c2{Xa,Xb,Xc)\" 


A pure valence quark description would only include the 
first term, without accounting for its Q^-dependence, 
which is governed by renormalization group equations 
coupling Cl to C 2 , C 2 to ci and C 3 , and so forth. Gen¬ 
erally speaking, for low values of Q^, i.e. < 1 GeV^, 
and most hadrons, with the possible exception of pions, 
the valence quark state will be the largest term in the 
expansion, but the contributions from other terms may 
not be negligible. 

This goes hand in hand with a change of the effective 
quark mass. While the assumption of quasi free partons 
with current quark masses less than 10 MeV for light fla¬ 
vors is valid for > 1... 2 GeV^, we know from lattice 
calculations and other considerations that the dynamical 
quark mass m{Q‘^) increases rapidly for < 1 GeV^ 
and ap pro aches a constituent mass of order 300 MeV 

dlEElii 

For the sudden recombination of medium constituents 
into a fast hadron, the relevant scale is determined by 


J 


the average thermal momenta in the medium. The pre¬ 
cise value of depends on the process that is being 
considered, but generally the scale falls into the range 
(ttT)^ < < (27rT)^, where T is the temperature of the 

medium [l5( . For a quark-gluon plasma near the point of 
hadronization {Tc « 170 MeV), this implies that the rel¬ 
evant scale Q is in the range 0.5 ... 1 GeV. It is therefore 
natural to assume that the correct degrees of freedom for 
recombination act like constituent quarks and that the re¬ 
combination probability is dominated by the lowest Fock 
states. 

Nevertheless, even for very massive constituents, 
higher Fock states should be present, even though we 
can not calculate their contribution from first principles. 
We want to discuss in the following, how an admixture 
of higher Fock states alters the recombination formalism. 

In the Boltzmann approximation, the probability for 
finding a quark (or gluon) with momentum k = xP and 
energy E\^ in the medium is given by: 

Wq{x) = {q{x)\p\q{x)) = (3) 

where p denotes the thermal density matrix and masses 
are assumed to be much smaller than Pc. Hence for a 
state with n quarks or other partons we have 

{q{xa)q{xb) ... IpI q{xa)q{xb).. .) 

= w,{xa)wgixb) ... = e-G^+--.+-)WT (4) 

The emission probability for the single Fock space com¬ 
ponents of a meson then is 


Wqq = / dXadXbS{Xa+Xb-l)\ci{Xa,Xb)\‘^ {q{xa)q{xb)\p\q{xa)q{xb)) =Cie (5) 

^0 

Wggg = / dxadxbdxcS{xa Xb X c - l)\c2{xa, Xb , X c)\^ {q{x a) q{xb)9 {xc) \p\ q{x a) q{xb)9{x c)) =^ 26 “^/^, ( 6 ) 
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etc. The delta functions expressing the condition that the 
sum of the light-cone momenta of all constituents add up 
to the light-cone momentum of the hadron ensure that 
the product of the Boltzmann factors of all constituents 
combine to the common factor containing only the 

hadron momentum. It is obvious that the same argument 
holds for any arbitrary complex Fock space component 
in the hadron wave function, as long as the medium is a 
thermal one. Note that we have suppressed quark flavor 
indices in the notation and factors arising from recombi¬ 
nation of color, spin, and flavor quantum numbers, which 
can be dealt with similarly but at the expense of a greatly 
complicated notation. 


Combining the contributions from all Fock space com¬ 
ponents, the probability for emission of a hadron with 
momentum P is given by: 

WiP) = Wgg + Wggg + Wgggg +■■■ = (7) 

where we applied the normalization condition ©• Our 
result shows that the probability of relativistic emission 
of a complex state by recombination from a thermal en¬ 
semble does not depend on the degree of complexity of 
the state. Per spin-flavor degree of freedom the emission 
of a baryon with momentum P is as likely as the emis¬ 
sion of a meson with the same momentum, as long as the 
particle masses are negligible compared with P = |P|. 
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We call this property the “egalitarian” nature of recom¬ 
bination from a thermal ensemble. 

Let us next explore how the elliptic flow of hadrons 
is affected by the presence of higher Fock states in their 
wavefunction. Here we assume that there are no space- 
momentum correlations affecting the calculation |l(il If7|| 
and we limit our discussion to sufficiently small values 
of the elliptic flow parameter V 2 , so that nonlinear cor¬ 
rections to the additivity rule of the flow of the con¬ 
stituents can be safely neglected |3|. The expression 
for the hadronic elliptic flow U 2 ^^(P) in terms of the 
quark/gluon elliptic flow V 2 {k) is: 



where v denotes the different Fock space configurations 
and a enumerates the constituents of the hadron in each 
configuration. Here it is important to be in a light cone 
frame where the masses of the particles can be neglected. 

In the hydrodynamic regime, where V 2 {k) = ak, one 
finds again trivially (P) = aP using @ . The elliptic 
flow of all hadrons then follows the same universal line. 
In the saturation regime, the relationship is more compli¬ 
cated and one may suspect that higher Fock states alter 
the V 2 of the hadron. 

For the derivation of the familiar scaling law for elliptic 
flow it is assumed that the wave function of the hadron 
is narrow: all partons in a Fock state carry roughly equal 
momentum Xi « l/n-i,, where rij/ is the number of par- 
tons. The scaling law follows when only the lowest Fock 
state with ni partons is taken into account 0 

= niV2iP/ni) (9) 

The experimental data is described very well by this 
equation HQ. We can easily generalize to higher 
Fock states in the limit of a very narrow wave function 
“ ^)\cu{xa)\^ « Cu Hi S{xi - l/n^). Then 

v‘''^\P) ~'^C„n„V2{P/n„) (10) 

I' 

and the scaling law is apparently violated by the contri¬ 
butions from higher Fock states. In principle, this vio¬ 
lation should be visible in a scaling analysis. The data 
are usually plotted with scaled axes Pt/tii and U 2 /ni, 
where ni = 2,3 is the valence quark number for the 
hadron. Equation m implies that the scaled elliptic 
flow for mesons and baryons, respectively, is 

(M) 

u 

(B) 

vf\p) = (SpM®)) . (11) 


Clearly, if all Q except the lowest Fock states, for 

which = 1 and = 1, the scaled elliptic flow 

curve is the same for mesons and baryons. This is what 
has been found in the data and has been interpreted as 
evidence that the scaled curve reflects the partonic ellip¬ 
tic flow before hadronization: u*^^)(p) = {;(®)(p) = V 2 (p). 

Can higher Fock states be present despite of the scaling 
law holding? We notice that for states with the same 
number s of sea quarks or gluons the ratio of prefactors 
riylni for mesons and baryons in Hill is 


3 s -1- 2 
2 s -1- 3 


( 12 ) 


and the ratio of prefactors inside the argument, multi¬ 
plying the momentum p is the inverse of that. This ratio 
is always between 1 for the pure valence state and 3/2 
for states with an infinite number of sea quarks and glu¬ 
ons. This is not a large variation, and there are two other 
effects that are also important. First, we note that the 
prefactors in front of V 2 and the momentum p in each 
line of cu tend to cancel if V 2 is a rising function. Sec¬ 
ond, with increasing number of sea quarks and gluons the 
average momentum of each parton decreases. For large 
Fock states we expect the partons to be in the hydrody¬ 
namic regime where V 2 is a linear function and scaling 
holds for arbitrary parton numbers. 

How large can the scaling violation be expected to be 
in practice? In order to present a numerical estimate, 
we consider the lowest Fock state and the next higher 
one, with an additional gluon, for mesons and baryons 
respectively. In focusing on the one-gluon admixture, we 
assume that the weights of higher Fock states, includ¬ 
ing two or more additional gluons or one or more quark- 
antiquark pairs, die out rapidly with growing complexity. 
For this analysis we abandon the simplistic (5-function 
shaped wave functions. We use the following more re¬ 
alistic light-cone distribution amplitudes for mesons as 
guidance: 


C^^\xa,Xb) 

y/XaXb 

{Xa,Xb,Xg) 


= 


'y(M) 


2V^’ 


^XaXbXg 

and for baryons: 
.(B) 


„(M) XaXbXg 
^ 12^10010’ 


{Xa^XbiXo) /^(B) XaXbXc 


y/XaXbXb 


{Xa,Xb,Xc,Xg) 

yfXaXbXcXg 


r^) XqXbXcX^g 

^ ^ 1680v^4^’ 


(13) 


(14) 


Here the variable Xg denotes the light-cone momentum 
fraction of the additional gluon in the hadron wave func¬ 
tion. The model wave functions on the right-hand side of 
m and da are standard forms of higher twist distribu¬ 
tions for hadrons M The additional kinematic factors 
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FIG. 1: Scaled elliptic flow of mesons and baryons (dotted 
lines) compared to the quark flow (solid line). The differ¬ 
ence between the quark and scaled hadron flow is due to the 
smearing by the internal hadron wavefunctions. 


FIG. 3: Scaled elliptic flow of mesons and baryons (dotted 
lines) compared to the parton flow (solid line), but now for 
hadron wavefunctions with a 30 percent component contain¬ 
ing an additional gluon. 
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FIG. 2: Scaled elliptic flow of baryons and mesons as calcu¬ 
lated from the quark flow compared to data for A-hyperons 
and kaons. 

on the left-hand side arise because of our non-standard 
probabilistic normalization of the states. 

The constants Ci give the relative normalization of the 
two Fock states. In the following, we shall consider two 
cases: (*) lowest Fock state only (Ci = 1,C2 = 0) and 
{ii) moderate \qqg) or contribution (Ci = 0.7, (72 = 
0.3). It is important to recognize that even in the first 
case, although valence quark scaling is well realized, the 
scaled elliptic flow function of the hadrons differs 
from the input function V 2 describing the quark flow. 
This is illustrated in Fig. Q] which shows the scaled el¬ 
liptic flow of mesons and baryons (dotted lines) together 


FIG. 4: Scaled elliptic flow of baryons and mesons as cal¬ 
culated from the parton flow, utilizing hadron wavefunctions 
with a 30 percent component containing an additional gluon, 
compared to data for A-hyperons and kaons. 

with the quark flow (solid line). The difference between 
the parton and scaled hadron flow is due to the smearing 
by the internal hadron wave functions. Figure |2] shows 
the scaled hadron elliptic flow in comparison with the 
scaled STAR data 0 for neutral kaons and A hyperons. 
The quark flow function was chosen so that the theoreti¬ 
cal values obtained by use of Eq. 0 fit the data. The fig¬ 
ure demonstrates that a single input curve for quarks can 
describe the elliptic flow of both, mesons and baryons. 

Figures 13 and 01 show the same quantities for the case 
{ii}, i. e. for hadron wave functions with a 30 percent 
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FIG. 5: Scaled elliptic flow of mesons and baryons (dotted 
lines) compared to the quark flow (solid line). The differ¬ 
ence between the quark and scaled hadron flow is due to the 
smearing by the internal hadron wavefunctions. 

component containing an additional gluon. We assume 
that the elliptic flow for the gluons is the same as that for 
light quarks. Note that the parton elliptic flow has been 
slightly adjusted to yield a good fit to the data again for 
this case. As anticipated, the scaling is not as good as for 
the lowest Fock state alone but, as Fig. |3 demonstrates, 
the scaling is only modestly violated. The largest viola¬ 
tion occurs for high transverse momenta. However, we 
also expect deviations to start in this region because frag¬ 
mentation starts to replace recombination as the domi¬ 
nate hadronization mechanism. The present precision of 
the data, as seen in Fig.01 is insufficient to observe the 
deviations from the scaling law and a gluon contribution 
of 30% is compatible with the data. 

Since the parton elliptic flow can be adjusted after in¬ 
cluding higher Fock states, the real test is the systematic 
deviation between baryons and mesons. In Figure [3 we 
show the relative difference ( 02 ^^— 02 ^^)/( 02 ^^be¬ 
tween the scaled meson and baryon elliptic flow for three 
different sizes of the higher Fock state component (0%, 
30%, 50%). In all cases, baryons have a slightly larger 
scaled V 2 than mesons at small momenta. This effect is 
likely to be overwhelmed by the influence of mass differ¬ 
ences, which have been neglected in the sudden recom¬ 
bination model. At larger momenta, the scaled meson 
V 2 is slightly larger. As stated, these differences are well 
within present experimental errors. 

In summary, we have investigated the effects of higher 
Fock states on the recombination of hadrons from ther- 
malized quark distributions. We showed that the yield 
of relativistic parton clusters is independent of the num¬ 
ber of partons in the cluster. Therefore, hadron spectra 
remain unaffected by the inclusion of higher Fock states. 


One important implication is that gluon degrees of free¬ 
dom could be accomodated during hadronization. They 
simply become part of the quark-gluon wave functions of 
the produced hadrons, but remain hidden constituents 
because the commonly produced hadrons do not contain 
valence gluons. 

On the other hand higher Fock states introduce devi¬ 
ations from the scaling law for elliptic flow. We showed 
that an additional 30% contribution from gluons is com¬ 
patible with the existing data on elliptic flow from RHIC. 
We emphasize that the interpretation that elliptic flow 
data from RHIC proves the existence of quark degrees of 
freedom in the bulk matter produced in the heavy ion 
collision is still valid. However, the connection of the 
measured elliptic flow to the quark elliptic flow might be 
less straight forward than anticipated. 
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